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We introduce a new non-Hermitian random matrix model for QCD with a baryon chemical poten- 
tial. This model is a direct chiral extension of a previously studied model that interpolates between 
the Wigner-Dyson and Ginibre ensembles. We present exact results for all eigenvalue correlations 
for any number of quark flavors using the orthogonal polynomial method. We also find that the 
parameters of the model can be scaled to remove the effects of the chemical potential from all ther- 
modynamic quantities until the finite density phase transition is reached. This makes the model 
and its extensions well suited for studying the phase diagram of QCD. 



Random matrix models have long been used as effec- 
tive models for a variety of complex quantum systems. 
More recently the class of non-Hermitian random matrix 
models has received much attention in systems including 
the fractional quantum hall effect Q], two dimensional 
charged plasmas |2j and quantum chaotic scattering (see 
for a review). Additionally one was introduced by 
Stephanov as a model for QCD with a baryon chemical 
potential |J]. Since then this model and related exten- 
sions have often been used to study the phase diagram 
of finite density QCD It is also useful for evaluating 
numerical algorithms that might be used to simulate fi- 
nite density QCD Q and in studying the density of the 
eigenvalues of the Dirac operator |jj . Despite the success 
of this model, exact results for the eigenvalue correla- 
tions have remained notably absent. Only recently has 
the eigenvalue density been obtained in the low energy 
(microscopic) limit by Splittorff and Verbaarschot ||. 

Another model for the eigenvalues of the QCD Dirac 
operator with a chemical potential was proposed by Ake- 
mann This model was constructed as a chiral exten- 
sion of a previously studied non-Hermitian matrix model. 
The main attraction of this model is that it can be solved 
exactly for all the eigenvalue correlations. However, the 
model is defined only in terms of the eigenvalues, and 
not from a matrix model based on the symmetries of the 
Dirac operator with a chemical potential. It was later 
shown that the partition function for this model agreed 
with that of Stephanov's in the microscopic limit |lflj . 
This agreement, however, only holds for a sufficiently 
small chemical potential as can be seen clearly by com- 
paring the microscopic spectral density @. 

Here we will consider a new matrix model for QCD 
with a baryon chemical potential which is also based on 
the form of the Dirac operator like in Stephanov's model, 
but can be solved for all spectral correlations like the 
model of Akemann. We will give exact results for the 
correlations with any number of quark flavors. We will 
also consider the results in the microscopic limit where 
the correlations should agree with QCD due to univer- 
sality. Lastly we look at the partition function and show 
that one can absorb the effects of the chemical potential 



by scaling the parameters of the model. This then pro- 
vides a physically accurate model that shows no change in 
thermodynamic quantities until the finite density phase 
transition is reached. 

The matrix model originally studied by Stephanov is 

with A a random complex matrix and /i is multiplied by 
the identity. This is based on the form of the Dirac oper- 
ator p + /i7o in a chiral basis. Here the matrix elements 
of 70 are chosen to be constant which is accurate when 
/i is very large, but it makes the model difficult to solve. 

Here we consider an alternate basis where the matrix 
elements of 70 are not taken as constant but are modeled 
with random matrices as 

Vn -\iA^ + ^ ) ■ (2) 

Here A and B are complex (N + v) x N matrices with v 
the topological charge. The QCD partition function with 
Nf quark flavors can now be modeled as 

. iff 
Z =j\f dAdB w G (A) w G (B) Y[det(V H + m 3 ) (3) 
J 3=1 

where Af is a normalization constant we will chose later. 
At low energies the results are universal and do not de- 
pend on the choice of the distribution for the matrix el- 
ements. We therefore make the simple choice of a Gaus- 
sian measure given by 

w G (A) = exp(-a NTtA^A) . (4) 

At fi = the parameter y/a = Y,V/2N is determined 
by the absolute value of the chiral condensate in the chi- 
ral limit (£) and the space-time volume (V). We will 
consider the /i dependence of a later. For now we will 
treat it as a constant and calculate the distribution of 
eigenvalues. 

In order to diagonalize the Dirac matrix we first make 
the substitutions C = iA + fxB and D = iA^ + fiB^ . The 
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partition function becomes 



becomes K v (x) ~ -\/7r/2xexp(— x) which gives 



r Nf 

Z=N' dCdD 102 (C, D) \\ det 



£> m., 



(5) 



with J\f' a dependent constant and the weight is 



w 2 (C, D) = cxp 



i ^ 



aN(l- n 2 ) 



Tr(CD + C ] D ] ) 



(6) 



The matrix I?// has exactly v zero eigenvalues and TV 
pairs of eigenvalues which we write as ±iz k - The non- 
zero eigenvalues of CD and DC are then equal to — z\. 
We now introduce the parameterization 



C 
D 



U{X + R)V 



(7) 



where U and V are unitary square matrices, X and R 
are complex (N + v) X N and Y and S complex N x 
(JV + v) matrices. The elements X. 
except when i 

i < j. The above decomposition is not unique and we 
must restrict the integration measure on U to the group 
V(N + v)/[V(l) N <g> U(v)] and on 1/ to U(7V)/U(l) Ar . 
The Jacobian for this transformation can be calculated 



and Yij are zero 
j and Rij and iSjj are zero except when 



i A (- 2 )i 2 n 



Xk 



(8) 



Zj) is the Vandermonde de- 



where A(z 2 ) = ]J i<jy 

terminant. Here we have used that XkVk = ~ z\ with 
and j/fc the diagonal elements of X and y respectively. 
With this parameterization the integrals over U, V, R 
and S become trivial and we are left with integrals over 
the complex numbers Xk and yu- We can change the in- 
tegral over yk to an integral over the eigenvalues Zk and 
then integrate out Xk to arrive at the form 



Z 



N f 



|A(, 2 )| 



l[w(z k )dz k Y[(zl + m 2 ) (9) 

fc=i j=l 



with N" another normalization constant and 

W (z) = \z\^e^( aN{l -/\ z* + z*>) 



4/i 2 



x K v 



aiV(l + /i 2 ) 
2/x 2 1 



(10) 



where K v is a modified Bessel function. 

The eigenvalue representation of the partition function 
© can be related to the model considered previously 
by Akemann. For large argument the Bessel function 



w(z) 



x exp 



aiV(l + M 2) l z l 
oJV(l + /i 2 ) 



2v+l 



(11) 



where we have set r = ( 1 — /x 2 ) / ( 1 + ^i 2 ) . Thisisthe same 
model considered in 9] . The two models agree when 



a7V(l + M 2 )kmm| 2 >2A< 2 



(12) 



The smallest eigenvalue z m i n has an average size around 
7r/(£V) at /i = 0. It is then convenient to consider the 
above condition in the microscopic |l ll| and weak non- 
Hermiticity |l2| limits given by taking N — > 00 while 
holding rjk = ^Vz k , — 2Nfi 2 and a fixed. The above 
condition then becomes 



\Vrmn\ 2 > 



(13) 



In the microscopic limit the two models agree only for 
sufficiently small fj, s . Since the new model is based di- 
rectly on the Dirac operator we expect that l(TU)) is the 
correct eigenvalue weight for QCD at low energies even 
for larger values of /i s . 

QUENCHED SPECTRUM 

Given the expression for the partition function in terms 
of the eigenvalues of the Dirac operator we can now solve 
for their correlations. This was done in for the model 
with weight (|lll) at Nf = using the method of orthog- 
onal polynomials extended to the complex plane. There 
it was shown that the corresponding orthogonal polyno- 
mials are the Laguerre polynomials. From that all spec- 
tral correlations can be determined analytically. It turns 
out that the orthogonal polynomials corresponding to the 
weight (|10fl are also Laguerre. We can therefore proceed 
in similar fashion. 

One can easily verify that the monic polynomials 



aN 



k\ LX 



( aNz 2 



are orthogonal with respect to the integral 



Pk( z ) Pi( z *) w ( z ) dz = h k hi 



Re(z)>0 



with 



hi = 



it fi {i + ii z y k+v k\ {k + v)\ 



2{aN) 



2k+v+2 



(14) 



(15) 



(16) 



Here we restricted the range of integration to the positive 
real half of the complex plane. This avoids double count- 
ing pairs of eigenvalues with opposite signs and ensures 
agreement with previously obtained results at /i = 0. 
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The 77-point spectral correlator can be written as 

p° N (z!, ...,z n )= det [K.° N (zi, Zj)] (17) 

l<i,j<n 



in terms of the kernel 



N-l 



K? N (x,y) = ^w{x)w{y) £ p° k (x) P ° k (y*) . (18) 



fe=0 



We can also evaluate the result in the microscopic limit. 
The microscopic limit of the 77-point correlator is 



P° s {vi,---,Vn) = Jim 



1 ( Vi 
-pN I 



Vn 



at^oo (£y) 2 ™ VEV'"' 



(19) 



with a kept fixed. This limit can be taken following the 
procedure in [i| . The result can be written as a determi- 
nant similar to (|17f) but with the microscopic kernel 



\xy 



ATrp 2 s (xy*y \ 



x e 



e~ 2 ^ ] J v (xVi)J v (y*Vi)dt 



o 



The spectral density is obtained by setting x = y in the 
above expression and agrees with the result given in || 
based on the model using Q). This agreement is expected 
since both models are based directly on the form of the 
Dirac operator. Only the microscopic results are univer- 
sal, though, and the macroscopic correlations of the two 
models do not necessarily agree. 

UNQUENCHED SPECTRUM 

We can also extend the results to an arbitrary number 
of flavors using an iterative procedure similar to the one 
used for the case of p = in |l3j. However since the 
mass term is only a function of z and not z* we can't 
use a single set of orthogonal polynomials. Instead we 
construct a set of biorthogonal polynomials for a flavors 
that satisfy the biorthogonality condition 

J p a k { Z ) q ?{ Z *)w(z)X[( Z 2 + m })dz = h a k 5 kl . (21) 



Re(z)>0 



3=1 



The monic polynomials can be defined iteratively by 
Pl+\^)pl- 1 {im a )-pl-\z) P a k - + \{im a ) 



p k {im a ) (z 2 + TO^) 
K' 1 ^ P^iima) _i/ *\ 



■J^g^OO (22) 



with the convention that <?°(z*) = p^(z*). The normal- 
ization coefficient can also be derived recursively by 



K 



Pk+l( im *) h a-l 
a— 1 / ■ \ k 

P k {tm a ) 



(23) 



An explicit formula for th e p olynomials for a general 
weight function appears in |14| . 

For convenience we will define the function 

a 

4>l{z)=p a k {z)X[{z 2 + rr%) . 

3 = 1 

The a flavor kernel can now be written as 



(24) 



K%(x, y) = lC a N (x, y) ^w{x)w{y) J] 



/ y 2 + m 2 - 



with 



N—l 



(25) 



(26) 



— n k 



fc=0 



We can get a recursion relation for this by substituting 
(221, © and (3J into (J2EI and rearranging the sums as 



iC a N (x,y)=iC a N - 1 (x,y) 



t a ^\im a ,y). (27) 



Repeated application of these relations gives an expres- 
sion for the a flavor kernel in terms of the zero flavor ker- 
nel and the functions p% +k (x) for < k < a — 1. Notice 
that the modified kernel JC^(x,y) vanishes at x = irrij 
for < j < a. It can therefore be written in the form of 
a determinant 



JC° N (x,y) p° N (x) ••• p° N+a -i( x ) 
KP N (imi,y) p° N {imi) ■■■ p° N + a -i( im i) 



iC%(im a ,y) p° N (im a ) ■■■ p° N+a _ 1 {im a ) 



(28) 



The denominator was determined by requiring that the 
coefficient of K? N (x,y) is unity. 

We can also take the microscopic limit of the above 
kernel. For this we also scale the masses with the inverse 
of the volume keeping 777^ = Y>Vm k fixed. This limit 
must be taken carefully and can be done with a procedure 
similar to that in |15j . The result is 



)C a s (x,y) = lC»(x,y) 



n 



K°(x,y) f = \^x 2 +m 2 



with tCg(x,y) equal to 



K°{x,y) J v {x) 



V J V +a-l(x) 



(29) 



K° s {irhi,y) J v {irhi) ■■■ (irhi) a J I/+0 _i(ij7ii) 



iC° s {im a ,y) Juiirha) ■■■ ( 



\(im k y l Ju+i-i(im k )\ 1<k4<a 



(30) 
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The function K° s {x, y) is denned as the integral that ap- 
pears in H20|) • This gives the general result for the micro- 
scopic spectral correlations of the QCD Dirac operator 
with a baryon chemical potential. 

PARTITION FUNCTION 

This model can also be used as an effective model for 
studying qualitative properties of the finite density QCD 
phase diagram. We can construct this model so that it 
has the correct physical behavior below the finite density 
phase transition. To see this we will look at the partition 
function which is simply given by 

N f N-l 

zZ'=/f'Nl]lm3l[hZ f . (31) 

3 = 1 k=Q 

Using the recursion relations l|23[l and (|22[1 the partition 
function can be written in terms of the quenched poly- 
nomials l|14|). The result is 

7 N S _ Kr N t - jj „ detigfc^jy, fe^iK)] 

3=1 V ' 

This result can also be obtained from the more general 
expression given in [lfij . 

Notice that the polynomials p° only depend on fi 
through the combination a/(l — /i 2 ). If we had started 
with a /i dependent a = (1 — /.i 2 ) ao and an appropriate 
normalization j\f then we could have made a partition 
function that doesn't depend on fj,. This means that all 
thermodynamic quantities would remain at their /i = 
values for fj, < 1. This agrees with the phenomenological 
expectation that the baryon number density should re- 
main at zero until a first order phase transition is reached. 
In this respect the new model is again an improvement 
over the previous one Q since that one shows an unphys- 
ical behavior in the baryon density as the chemical poten- 
tial is increased. Another matrix model with improved 
thermodynamic properties was studied in [l7j . however, 
that model has not been solved for the eigenvalues. 

At fj, > 1 the above scaling is no longer valid since the 
Gaussian integrals in the partition function would not 
converge. However, one can show that for /i > 1 and 
a > the chiral condensate in the unquenched theory 
vanishes in the large N limit. The model therefore does 
have a first order phase transition at /.i = 1 . Of course for 
large fi the Dirac operator should approach the determin- 
istic form /Lt7o given in (I). Thus this new model cannot 
be expected to correctly describe the eigenvalue correla- 
tions in the phase where chiral symmetry is restored. It 
is currently then a phenomenological model for QCD at 
fi< n c (here \i c = 1). 

Given this models success it would be interesting to 
include temperature into the model and study the full 



phase diagram as has done with previous matrix mod- 
els |5j- A more complete analysis of the phase diagram 
will be saved for a later publication. Another interesting 
extension would be to make both A and B be random 
banded matrices with a power law decay like what was 
considered for \i — in 18]. This was found to provide 
good agreement with an instanton liquid model for eigen- 
value correlations at larger energies (beyond the Thouless 
energy). This also may provide an accurate way to study 
the finite temperature chiral restoration transition. 

In conclusion we have introduced a new random matrix 
model for QCD with a baryon chemical potential that is 
based on the form of the Dirac operator in a chiral ba- 
sis and is exactly solvable for all eigenvalue correlations. 
We have presented the correlations for any number of 
quark flavors and given results in the microscopic limit 
where they should agree with QCD due to universality. 
This new model also has a physically accurate thermo- 
dynamic behavior for /x < \i c and could provide a basis 
for a phenomenological model of the finite density QCD 
phase diagram. 
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